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ABSTRACT 

We present spectroscopic observations of the quiescent black hole binary A0620-00 
with the the 6.5-m Magellan Clay telescope at Las Campanas Observatory. We mea- 
sure absorption-line radial velocities of the secondary and make the most precise de- 
termination to date {K2 = 435.4 ± 0.5 km s""'^). By fitting the rotational broadening 
of the secondary, we refine the mass ratio to g = 0.060 ±0.004; these results, combined 
with the orbital period, imply a minimum mass for the compact object of 3.10 ± 0.04 
Mq. Although quiescence implies little accretion activity, we find that the disc con- 
tributes 56 ± 7 per cent of the light in B and V, and is subject to significant flickering. 
Doppler maps of the Balmer lines reveal bright emission from the gas stream-disc 
impact point and unusual crescent-shaped features. We also find that the disc centre 
of symmetry does not coincide with the predicted black hole velocity. By comparison 
with SPH simulations, we identify this source with an eccentric disc. With high S/N, 
we pursue modulation tomography of Ha and find that the aforementioned bright 
regions are strongly modulated at the orbital period. We interpret this modulation in 
the context of disc precession, and discuss cases for the accretion disc evolution. 



Key words: 

close 



accretion, accretion discs — stars; individual: A0620-00 — binaries: 



o 



X 



1 INTRODUCTION 

A0620-00 (V616 Mon) is the prototype Soft X-ray Transient, 
a class of low-mass binary stars which exhibit infrequent but 
intense X-ray bursts (Gelino, Harrison, and Orosz, 2001). In 
1975 it became the bri ghtest X-ray nova ever detected, at ap- 
proximately 50 Crab (|Elvis et al.lll975l ). and it was the first 
nova to be identified with a black hole primary (McClintock 
& Remillard, 1986; hereafter MR86). MR86 measured an or- 
bital period of 7.75 hr and a radial velocity semiamplitude 
for the K-type secondary of 457 km s"'^, leading to a mass 
function f{M) = 3.18 Mq; estimates of K2 and f{M) have 
decreased slightly since then (i.e. 433 km s"^ and 3.09 Mq) 
(Marsh, Robinson, & Wood 1994, hereafter MRW94). Given 
this minimum mass, it is likely that A0620-00 is a black hole. 

A substantial amount of work has gone into the anal- 
ysis of A0620-00 in the last twenty years, with particular 
emphasis on ellipsoidal variations in the light curve and the 
contamination of the K-star flux by light from the accretion 
disc. As yet, no real consensus has been reached, mostly due 
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to the complexity of the lightcurves. While ellipsoidal vari- 
ations are obvious, they are highly asymmetric (Leibowitz, 
Hemar, & Orio 1998); the origin of t he asymmetry is un - 
determined. Modelling this lightcurve, iGelino et al.l (j200ll ) 
determined in inclination of 41±3°, invoking starspots to ex- 
plain the asymmetries. Shahbaz, Naylor, and Charles (1994) 
found a 90 per cent confldence interval of i =30-45° given 
the mass ratio of A0620-00, modelling their asymmetries 
with the bright spot where the accretion stream hits the 
disc. 

Lightcurve modelling is also complicated by the vari- 
ability of the disc itself. To quantify ellipsoidal variations, 
most authors assume the disc to be constant, and justify 
the claim by noting that A0620-00 is quiescent. They do not 
mention that estimates of the dis c contamination range from 
<3 per cent jGelino et aLlbOOll ) to < 50 per cent (MR86). 
The contribution from this disc is not only unclear, but ap- 
parently not constant. More than half of A0620-00's 58-year 
burst cycle has passed, and it is important to note that qui- 
escent does not mean inactive. We will argue that the vari- 
ability of the accretion disc cannot be neglected. In order 
to determine definitively the mass of the compact object, 
it is very important to understand the structure and varia- 
tion of the accretion disc. MRW94 made enormous progress 
towards this goal. In 2004, Shahbaz et al. (hereafter S04) 
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Table 1. Observation log 



Source 


Date 


# of Spectra 




A0620-00 


2006 Dec 14 


11 


873 


A0620-00 


2006 Dec 15 


23 


698 


A0620-00 


2006 Dec 16 


14 


433 


HILT 600 


2006 Dec 14 


1 


20 


HILT 600 


2006 Dec 15 


3 


20 


HILT 600 


2006 Dec 16 


3 


20 


HD 18298 


2006 Dec 15 


3 


5 


HD 7142 


2006 Dec 15 


5 


5 



Tcxp is the average exposure time for the source. 

noticed signatures of an eccentric disc not seen in previous 
Doppler maps, but lacked the phase coverage to verify their 
hypothesis. 

Therefore, as follow-up to the work of MRW94 and S04, 
and as part of a Doppler imaging survey of black hole and 
neutron star binaries, we undertook phase-resolved optical 
spectroscopy of A0620-00. In §2 we describe our observa- 
tional methods and data reduction. In §3 we measure the 
radial velocity of the secondary star, the system mass ratio, 
and attempt to quantify flickering. In §4 we present Doppler 
images of the accretion disc at several wavelengths, inves- 
tigate evidence for an eccentric disc, and report results of 
modulation tomography of the Ha line. We discuss conclu- 
sions from the variability of the disc and our Doppler maps 
in §5. 



2 OBSERVATIONS 

We observed A0620-00 with the Low-Dispersion Survey 
Spectrograph (LDSS3) at the f/4 focus of the 6.5-m Clay 
telescope at Las Campanas observatory on 2006 December 
14-16. We acquired 48 spectra using the VPH Blue grism 
and a long 0.75 arcsec slit. By shifting the slit 4° redward, 
we were able to observe Ha with the superior resolution of 
the Blue grism (2.3 A = 130 km s"^), covering 4250-7035 
A. To minimize the effects of atmospheric dispersion, we ob- 
served at parallactic angle. Our exposure times ranged from 
420 s to 1200 s, with an average of 660 s (a total of 8.81 
hours on the source). 

Each night we observed the flux standard HILT 600 
with the same instrumental setup as A0620-00. On all nights 
the seeing was generally comparable to our slit width, but 
occasionally spiked as high as 1.8 arcsec due to wind, and slit 
losses prevented precise flux calibration. On 2006 Dec 15 we 
also observed two K3/K4 stars, HD 18298 and HD 7142, as 
velocity standards, again with the same optica l setup. As the 
second ary is constrained to be later than K3V lFroning et al] 
we expect more accurate results for HD 7142, which 
is listed as K3/K4III in the SIMBAD database (HD 18298 
is listed as K3IIICN). For each pointing, we obtained com- 
parison HeNeAr arc lamp spectra after every 3-7 spectra, 
depending on the current exposure time. We list the obser- 
vations in Table 1. 

We used standard IRAF routines for basic data reduc- 
tion (zero-subtraction, flat-fielding, and spectral extraction). 
We extracted our spectra in multispec format, attempting 
apaU's optimal extraction using nominal LDSS3 gain, read- 
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Figure 1. The average normalized spectrum of A0620-00. The 
strongest features are, in order of increasing wavelength, H7, H/3, 
He I 5875, Ha, Hen 6678, and telluric absorption lines. We also 
show the scaled spectrum of HD 7142 for reference. 



out noise, and full-well values. The routine also performs 
standard extraction and generates error bars. A CCD de- 
fect running across our spectra prevented reliable fits to the 
spatial profile for optimal extraction, so we used the nor- 
mally extracted spectra instead. This choice did not signifi- 
cantly degrade our S/N, and should have a negligible effect 
on our presented results. After wavelength calibration, we 
passed the spectra and their errors to the software package 
MOLLY for cleaning and analysis. 



3 ANALYSIS 

3.1 The spectrum of A0620-00 

In Figure 1 we present the average normalized spectrum of 
A0620-00. The spectrum shows a number of strong features 
(originating from the disc) and many relatively weak K-star 
absorption lines. We overplot the K-dwarf template HD 7142 
for reference. We have extremely high S/N Ha and H/3 lines, 
both observed by MRW94, as well as lower S/N H7 and 
He I 6678A. The non-detections of Hell lines indicates the 
scarcity of ionizing radiation. Our Ha line is stronger relative 
to the continuum than in 1994 by approximately 50 per 
cent. This relative brightening is not surprising given the 
decreased fraction of light contributed by the secondary. We 
show a close-up of Ha in panel a of Figure 2. The profile 
has two strong symmetric peaks. It is also interesting that 
our H/3 line also shows two strong peaks (as opposed to a 
single-peaked line seen by MRW94). 

The other feature of note is located on top of the in- 
terstellar sodium doublet near 5890 A. We show a close-up 
of this line in panel b of Figure 2. The line suffers signifi- 
cant extinction by interstellar sodium, so it is not possible 
to identify its peak unequivocally, but it seems to be He I 
5875. It appears to have some structure, and trailed spectra 
suggest a double-peaked profile, but our attempts to correct 
for interstellar absorption and create a Doppler map (§4) 
were unsuccessful. 
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Figure 2. Close-ups of a) the H« line and b) the feature near 
5890A. Intensities are as in Figure 1. 



3.2 The radial velocity of the secondary star 

Following MRW94, we first measure the radial velocity of 
V616 Mon in order to estimate the mass function. Mask- 
ing out emission lines, telluric absorption, and interstellar 
sodium lines in our velocity standards, and normalizing, 
we cross-correlated our spectra against these K-dwarf tem- 
plates. The results of the following sections are presented in 
Table 2. 

As discussed in MRW94, this process requires several 
adjustments for A0620-00: rotational broadening and orbital 
smearing. The following analysis was completed for each 
template. We performed a preliminary broadening of the 
template with MRW94's value of i;sini = 83 km s"^. Then, 
after cross-correlating, we shifted our spectra into the rest 
frame of the template, which is essentially the rest frame 
of the secondary. We made 48 copies of the template and 
smeared each according to the formula s = 2'KViTi/P, where 
P is the orbital period and Vi and Ti are the radial veloc- 
ity and exposure time for the ith A0620-00 spectrum. With 
exposure times of up to 1200 s, it is necessary to correct for 
smearing because s is comparable to v sin i. We averaged the 
smeared templates, rotationally broadened the result with 
a value between 30 and 135 km s~^, and optimally sub- 
tracted the final template from the average object spectrum. 
We assume a limb-darkening coefficient e of 0.65 (Wade & 
Rucinski 1985, and references therein), but perform these 
calculations for e — 0.45 — 0.85 to evaluate our systematics. 
The optimal subtraction routine returns the fraction / of 
light contributed by the secondary and a value. We fit 
fourth-order polynomials to fo identify the appropriate 
/ and vsini. Results are included in Table 2. 

The measured rotational broadening of V616 Mon (80± 
2 for HD 18298 and 83 ± 2 for HD 7142) is consistent with 
MRW94, but / is not, even though / and v sin i were strongly 
correlated (i.e. we were able to identify them with the same 
X^ minimum) . MRW94 found that the secondary contributes 
~ 94 per cent of the light near Ha and ~ 85 per cent of the 
light near H/3. We observed both wavelengths simultane- 
ously, and find that the secondary contributes only 44 per 
cent of the light in B and V (near 5500 A). This conclusion 



is effectively independent of the limb-darkening coefficient. 
We will discuss systematic uncertainties on / in §3.4. 

Finally, we performed another cross-correlation of our 
spectra against the template, this time with the appropriate 
value of usini. The results are shown in Figure 3. Although 
we took arc exposures frequently to minimize the effects 
of flexure, we found arc scales to drift by ~0.5 A over the 
course of a few hours. As our dispersion (0.69 A pixel"^ 
) corresponds to 37.9 km s~^ pixel"^, it was necessary to 
compensate for arc drift. We did this by cross-correlating 
the telluric lines near 6900 A and shifting out the result- 
ing velocities (generally around 10 km s~^). In this way we 
guaranteed a common (heliocentric) rest frame for all 48 
spectra. 

We fit a function of the form V ^ K2 sm{2Tv(t-to)/ P) + 
7 for each template. 7 represents the systemic velocity in 
the template frame. Although we have very high S/N, the 
baseline (~7 orbits) was insufficient to make a reliable in- 
dependent determination of the orbital period, so we fixed 
P according to the ephemeris of MR86. To get an accu- 
rate measure of our statistical uncertainties, we performed 
10,000 Monte Carlo simulations for each fit, assuming each 
radial velocity to be distributed normally around the mea- 
sured value, increasing the errors from cross-correlation by 
a factor of three. For our final values, we averaged the re- 
sults from both templates and took the standard deviation 
of the combined distribution as the uncertainty. The result- 
ing uncertainty is not purely statistical, as the use of multi- 
ple templates includes some systematic errors (e.g. template 
mismatch). These results can be found in Table 2. 

All fits are comparable in quality and have excellent x^ ■ 
We find consistent values for K2 and excellent agreement in 
to ~ we observed inferior conjunction on the second night 
- and our measurements are independent of limb-darkening 
(e=0.65 is marginally preferred). We will consider the sys- 
temic velocities only briefly. The uncertainties in 7 are dom- 
inated by the systematic uncertainty (one flfth of a pixel, 
or ~7.6 km s~^). We were unable to find a cataloged ra - 
dial velocity for HD 18298, but from lMalaroda et al.l (j200lh . 
the radial velocity of HD 7142 is 32.8 km s"^ Given that 
MRW94 report a systemic velocity of 22 km s~^, we con- 
sider these fits to be accurate; the choice of template does 
not significantly affect K2, to, or wsini. In summary, we 
adopt vsini = 82 ± 2 and K2 = 435.4 ± 0.5 km s"\ In ad- 
dition, we find that inferior conjunction of the mass donor 
star, which defines the zero point to of our ephemeris, oc- 
curs at HJD (UTC) 2454084.69485±0.00005. This latest to 
is consistent with the ephemeris of MR86 (within 0.5 a). 



3.3 The mass ratio 

To measure the mass ratio, we use Paczynski's (1971) for- 
mula relating the rotational broadening and K2 to q for 
Roche lobe-filling stars. 



= 0.462[(l + g)"g] 



1/3 



(1) 



K2 

This gives q — 0.060 ± 0.004, where the uncertainty in- 
cludes variation between templates. This value is consis- 
tent with MRW94's measurement (by the same method) of 
0.064 ± 0.01, and with their value of q = 0.067 ± 0.01, ob- 
tained by calculating models on a grid over the Roche lobe. 
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Figure 3. Radial velocity of V616 Mori. Two cycles are shown for clarity. Filled circles correspond to the template HD 18298, and open 
circles to HD 7142. Note from Table 2 that the main difference is in 7. 



Table 2. Fits to rotational broadening and radial velocities. 



Template e vsini f K2 (km s"^) to - 2454084 7 (km s"!) q /(Afi) (Mq) 

(km s-i) ±2 ±0.001 ±0.4 ±4E-5 ±8 ±0.004 ±0.03 



HD 18298 


0.45 


77 


0.441 


435.8 


0.69487 


11 


0.051 


3.06 


1.06 




0.55 


78 


0.441 


435.8 


0.69487 


11 


0.052 


3.07 


1.06 




0.65 


80 


0.440 


435.8 


0.69487 


11 


0.057 


3.09 


1.05 




0.75 


80 


0.441 


435.8 


0.69487 


11 


0.056 


3.09 


1.06 




0.85 


81 


0.441 


435.8 


0.69487 


11 


0.058 


3.10 


1.06 


HD 7142 


0.45 


80 


0.433 


435.0 


0.69484 


-9 


0.056 


3.07 


1.03 




0.55 


80 


0.433 


435.0 


0.69484 


-9 


0.057 


3.08 


1.03 




0.65 


83 


0.449 


434.9 


0.69484 


-9 


0.063 


3.11 


1.02 




0.75 


83 


0.433 


435.0 


0.69484 


-9 


0.061 


3.10 


1.03 




0.85 


84 


0.433 


435.0 


0.69484 


-9 


0.064 


3.12 


1.03 



including gravity darkening and quadratic limb darkening. 
Although it might be possible with our improved resolution 
to distinguish between the Paczynski approximation and the 
grid models, since the grid model correction was far from 
significant at their la level, we opt to take our results as 
accurate. This choice is validated by the weak (at best) de- 
pendence of our measurements on limb darkening, and we 
report uncertainties large enough to account for any system- 
atics in e. 

Given the orbital period, K2, and q, the minimum 
masses for both objects are: 



where the true mass goes like Mmin/ sin'^ i. As noted in many 
papers (MRW94; Shahbaz et al. 1994), the mass of a maxi- 
mally rotating neutron star with the stiffest equation of state 
is 3.2 Mq; if causality is the only constraint, t he absolute 
upper limit is 3.76 Mq (iFriedman fc Ipse i] 19871 ) . Using the 
constraints 39° ^ i 75° (|Gelino et all 120011 )7 we find 3.4 



Mq < Ml ^ 12.6 Mq, with the most likely value at 11.1 
Mq. We thus improve the precision of the mass function, 
but it remains a possibility that A0620-00 is a stiff, massive 
neutron star. 



f{Mi) = 3.10 ± 0.04 Mq 
/(Ma) = 0.19 ± 0.02 Mq, 
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Figure 4. a) Equivalent width of the H« line, b) Secondary light 
fraction as measured from 5000-5800 A and 6000-6500 A for the 
template HD 7142. Differences from HD 18298 are less than la. 
These plots indicate a flickering source of light. The solid, dotted, 
and dashed lines are the predicted light fraction for a constant 
disc and ellipsoidal modulations, for inclinations of 40°, 60°, and 
80°, respectively. See Equation 2. 



3.4 Ellipsoidal variations and the variability of 
the disc 

Part of the difficulty in measuring the parameters of A0620- 
00, and a possible source of systematic error, is the tidal 
distortion of the secondary. We assume, for example, in our 
radial velocity measurements, that the secondary's centre 
of light coincides with its centre of mass. For a star filling 
its Roche lobe with some gravitational and limb dar kening, 
this a ssumption is not likely to be valid. Following ISternel 
l|l94ll ). we checked for this distortion in our data by fitting 
an extra — JsTeii sm(An(j) ) term to our radia l velocity curve. 
Orbital elements from iGelino et al.l (|200ll ') lead to a pre- 
dicted i^cii = 3.0 km s~^, but we found Kdi ~ 0.7 km s^^ 
(less than the error bars on each radial velocity), and an 
insignificant improvement to . It is possible that our data 
have insufficient S/N to make this measurement. This par- 
ticular effect is small, but we must be careful not to dismiss 
the systematic uncertainties from ellipsoidal effects. 

Ellipsoidal variations, particularly in the lightcurve, re- 
ceive a great deal of attention because of the constraints 
they place on the orbital parameters, especially the inclina- 
tion. In practice, one usually attributes all variation to the 
secondary (assuming the disc to be constant). However, the 
fact that the fraction of light contributed by the disc changes 
by a factor of nine over ten years casts significant doubt on 
that particular assumption. In turn, this evolution must be 
accounted for when interpreting lightcurves spanning a long 
period of time. 

Given that we have good spectral resolution and rela- 
tively good phase coverage, we have attempted to charac- 
terize the ellipsoidal variations of the source during our ob- 
servations. While it is most common to address variations 
in the lightcurve, effects of tidal distortion should also be 
apparent in the secondary light fraction and the line equiv- 
alent widths. In the rest of this section we quantify these 
variations. 

If the disc is a constant diluting source of light, then con- 



tinuum variations, and therefore variations in the equivalent 
widths of lines, maybe be attributed to the tidal distortion 
of the secondary. In Figure 4a we show the equivalent width 
of the Ha line as a function of orbit al phase. The asymm et- 
ric ellipsoidal variations observed bv lCelino et al.l (|200lh are 
obvious, along with noticeable scatter. Given a S/N of ~ 95 
at Ha, it seems reasonable to interpret this as fiickering, al- 
beit undersampled, rather than noise. We will discuss the 
equivalent widths in more detail shortly. 

Figure 4b shows the secondary light fraction / as a func- 
tion of orbital phase. We can evaluate the variability of the 
disc with a back-of-the-envelope calculation of ,f{(j)). We as- 
sume the secondary continues to exhibit ellipsoidal varia- 
tions of 0.077 mag (MRW94 using i = 40°), or flux varia- 
tions of 7 per cent, and that the disc is constant in time. 
Since the disc contributes 56 per cent of the light, it should 
be roughly 22 per cent brighter than the secondary. Then 
in units of the mean secondary flux, F2 = 1 — O.O7cos(47r0) 
and Fx = 1.22, and 



/ = 



1 - 0.07cos(47r<;/>) 
2.22 - 0.07cos(47r(/> 



(2) 



The inclination-dependent amplitude of ellipsoidal varia- 
tions is tabulated in MRW94. Equation 2 is plotted along 
with / for inclinations of 40, 60, and 80°. We have matched 
the mean relatively well, and to some extent the functional 
form, but it is obvious that our assumptions do not hold. 
First, the disc appears anomalously bright between =0.7- 
0.9. As this phase interval was observed on the second night 
only, we cannot speculate if the dip in / indicates a flare 
or a long-lived bright region of the disc, like a warp, visible 
only at this phase. 

The second and more troubling discrepancy is the large 
amplitude of variation of /. As is apparent from Figure 4b, 
excluding (f> =0.7-0.9, higher inclinations are preferred, even 
those which are ruled out by the lack of eclipses in this 
so urce. Accepting mom entarily the inclination determined 
by iGelino et al.l ( 200 ih , this figure illustrates very clearly 
the significance of the assumption that the disc is a con- 
stant source of light; by requiring a constant disc, we could 
overestimate the inclination by several tens of degrees. Now 
it is reasonable to assume that the secondary has not evolved 
substantially in the last twenty years, so we must conclude 
that the discrepancy between the predicted and observed 
light fractions is related to activity in the disc. For example, 
a component of disc light modulating at the orbital period 
could reproduce the effect easily. 

Unfortunately, there is no a priori way to determine 
in advance the viability of an inclination measurement. 
SMARTS data from the last ten years show that A0620-00 
goes through periods of quiescence, where ellipsoidal modu- 
lations are observed, and periods of erratic variability which 
gradually swamps the smooth component (Cantrell & Bai- 
lyn 2007, private communication). In these active periods, 
the source is extremely variable on timescales from minutes 
(our observations) to years (SMARTS), and the variability 
is highly phase-dependent, so that ellipsoidal modulations 
cannot be reliably measured; the current active state has 
persisted since December of 2003. Associating this variabil- 
ity with the disc, we warn against measurements of i during 
such periods. We can expect, furthermore, a strong corre- 
lation between flickering and the light fraction. There is no 
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reason / cannot be small and constant simultaneously, but 
an increase in light from the disc can only mean an increase 
in accretion activity, for which flickering is highly probable. 

This is admittedly a grim assessment of the situation, 
but lightcurve estimates of the inclination made during pe- 
riods of variability or non-negligible contamination by disc 
light are unreliable. The variation in / itself is a final source 
of systematic uncertainty. We have a very precise measure- 
ment of the dilution fraction for the average spectrum, but 
that precision is only meaningful if / is not variable. Instead, 
we take the standard deviation in f{<j)) as our uncertainty: 
/ = 44 ± 7 per cent, and summarize the discussion above by 
reminding the reader that quiescence is not inactivity, and 
ought not be used to justify invalid assumptions. 

Consider now the equivalent width. Since the contin- 
uum varies slowly over our lines, it is trivial to show that 
the line equivalent width is given by 



EW (A) 



(3) 



where _EWdisc is the ratio of flux integrated over the line to 
the disc continuum, and / is the secondary light fraction. If 
the disc varies uniformly or not at all, EWdisc should be a 
constant for any given line. We can then attribute variations 
in the equivalent width to variations in /, which would be 
ellipsoidal in nature. However, measurements of the equiva- 
lent width with the secondary subtracted should reveal the 
reliability of this shaky assumption. 

The equivalent width of our unsubtracted Ha line (open 
circles) is shown in Figure 5; it should be compared to the 
same plot from MRW94. To facilitate this comparison, we 
have overlaid their best-fitting line. We find good agreement 
with their phasing and amplitude, but an increase in average 
equivalent width of approximately 13 A. As noted earlier, 
the H-alpha line is brighter relative to the continuum than 
it was in 1994, so this increase is reasonable. Also shown in 
Figure 5 (filled circles) is the equivalent width of the same 
line after subtraction of the secondary. To achieve this re- 
sult, which indicates the magnitude of fluctuations in the 
disc, we normalized, broadened, and smeared the template, 
subtracted /(</>) times the template from each spectrum, and 
set the continuum to one. The dip between phases 0.7 and 
0.9 corresponds to the possible flare seen in the light frac- 
tion. 

We also note that the subtracted modulations are in 
phase with the unsubtracted equivalent widths. As the sec- 
ondary is expected to contribute more light at longer wave- 
lengths, the phasing may be an artefact of improper sec- 
ondary subtraction, because the strongest contributions to 
/ come from 5000-5800 A. However, this effect should be 
small because absorption lines up to 6500 A were used in the 
measurement. Even if this is the case, the scatter in EWdisc 
cannot be explained by the combined noise in / and EW. 
We suggest that the most probable extra source of scatter 
is physical variability of the disc, i.e. flickering. In this case, 
the fact that the subtracted equivalent width is not con- 
stant suggests that the disc fluctuations cannot be spatially 
uniform. 

Our measurements of flickering are confirmed by nearly 
simultaneous observations with all four IRAC bands on 
Spitzer and the 1.2m FLWO telescope (McClintock 2007, 
private communication). These observations, taken in the 



■t 100 




Figure 5. Measured equivalent widths for the Ha line from 
A0620-00. Open circles correspond to unsubtracted spectra; filled 
circles were calculated after subtracting the template HD 7142. 
Again, differences between templates were insignificant, and the 
errors are smaller than the circles because of the extreme signal- 
to-noise. We have plotted MRW94's best fit for comparison. 



last week of November 2006, show strong erratic variability 
which is well-correlated between telescopes. Neither ellip- 
soidal variations nor the orbital period is obvious through 
these fluctuations. The conclusion is clear: the assumption 
that the disc is constant is not valid, though it may be a 
good approximation if the disc contrib ution to the c ontin - 
uum is negligible (for example, when iGelino et al.l (|200lh 
found / > 97 per cent and i = 41 ± 3°). To investigate 
further the spatial and temporal intensity of the disc, we 
present the results of Doppler tomography and modulation 
tomography in § 4. 



4 DOPPLER IMAGING 

In this section we discuss the results of Doppler tomogra- 
phy of the emission lines from A0620-00. To create each 
image, we subtracted the scaled HD 7142 spectrum, per- 
formed a linear fit to the continuum around each line, nor- 
malized, and set the surrounding continuum to zero (re- 
quired by DOPPLER). MRW94 express some uncertainty as 
to the propriety of interpolating over the Ha line in the 
template spectrum. Doing so removes an image of the donor 
star in the map, and it is not clear which is the appropri- 
ate choice. However, we found that Doppler maps including 
the light from the secondary do not show donor emission, 
so we chose to interpolate over the line. Then we binned the 
spectra to a uniform velocity scale and passed them to the 
DOPPLER routine, which computes the maximum entropy 
image that reproduces the observed line profiles during the 
course of an orbital period . For an excellent summary of the 
method, see lMarshI l|200lh . 



4.1 The maps 

In Figures 6-9 we present Doppler tomograms of the strong 
emission lines in A0620-00. In each figure, the left column 
is, from top to bottom: the observed trailed spectrum, the 
Doppler map, and the fitted data. The right column is the 
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observed data minus a simulated symmetric part, the asym- 
metric part of the Doppler map, and its fitted data. In the 
maps, we plot the Roche lobe of the secondary, the ballistic 
trajectory of the gas stream (lower line), and the Keplerian 
velocity of the disc along the stream (upper line) . 

In Figure 6 we show the maps of Ha. The most obvi- 
ous feature is the bright spot, which corresponds to the gas 
stream impact point. As found by MRW94 and S04, in our 
maps the gas stream trajectory and the Keplerian disc ve- 
locity along t he stream (the two lines plotted in the map 
panels) bracket the bright spot. We can interpret this as 
post-shock emiss ion, originating so mewhat inside the outer 
edge of the disc (|Marsh et al.lll990l ). Using the ballistic tra- 
jectory, we locate the spot at r = 0.6 ± 0.057?li. If the disc 
velocities are Keplerian and the spot moves with the disc, 
we find the outer edge near r — 0.45±0.05iiLi. This is some 
indication of our systematic uncertainties, but we shall sug- 
gest shortly that the larger disc is more likely. 

The other features of note are the two crescents at ~ 7 
o'clock and ~ 2 o'clock. We shall address their origin shortly. 
The fitted data are rather messy, mainly due to the pres- 
ence of those features and the strong flickering. Because the 
crescents and the bright spot so thoroughly dominate the 
image, the symmetric part is too bright, and we oversub- 
tract to make the asymmetric part. Still, we do detect some 
emission from the bright spot. The trail itself is not partic- 
ularly sinusoidal: it is more of a zigzag than an S-wave. The 
discrepancy is best seen in the simulated asymmetric trail 
(bottom right, Figure 6), but is apparent in the top panels 
as well. In the data panel (top left), the trail is diagonal 
between (j> = —3 and (j) — —2.5, but nearly horizontal at 
<j) = —2. It is difficult to interpret this as the S-wave of a 
circular orbit. If the orbit is not circular, it may be a com- 
bination of the ballistic trajectory of the stream and the 
motion of the disc, or it may be an elliptical orbit. 

In Figure 7 we show maps of H/?. Again, we see the cres- 
cent features; their shape is more apparent despite the lower 
S/N at H/3. For the same reason, we have less trouble com- 
puting the symmetric part, and our asymmetric part shows 
some emission along the stream. While the corresponding 
trail does look nicer than its Ha counterpart, the two seem 
to be consistent. In Figure 8 we show maps of H7, and in 
Figure 9 we show maps of He I 6678. These two maps, at 
much lower S/N, reveal bright spot emission with some con- 
tribution from the stream itself. 

4.2 The eccentric disc of A0620-00 

In the previous sections, we have presented evidence for a 
bright flickering disc which contributes more than half of 
the light from the source in the optical. It is already clear 
that the "no intrinsic variability" assumption of Doppler to- 
mography is not satisfied by A0620-00. We should therefore 
exercise caution in interpreting our Doppler maps. Tomog- 
raphy is a robust technique, but the physical relevance of its 
results depends on the extent to which its assumptions are 
violated. For example, we assume that all velocity vectors 
corotate with the binary. However, if the disc is large enough, 
it may reach an orbital resonance with the secondary and 
become eccentric by tidal distortions. It will then precess, 
even in the corotating frame. For the mass ratio of A0620-00 , 
the dominant resonance is 3: 1 (|Whitehurst fc Kind[l99ll) . 



Recent SPH simulations bv lFoulkes et al.l \2004 ) show that 
the manifestations of eccentric discs are bright emission be- 
tween the gas stream and the Keplerian disc velocity along 
the stream, non-sinusoidal S-waves, crescent-shaped features 
in Doppler maps, and the shifting of the map centre of sym- 
metry away from (0,--R'i). 

S04 observed the first three features and invoked an ec- 
centric disc, but they also point out one possible objection: 
MRW94 found a clean circular disc. However, we have expo- 
sure times comparable to MRW94 with a larger telescope; 
our higher signal to noise may have enabled us to detect 
this phenomenon. The other explanation is that the cres- 
cent was absent in 1994, and the source has changed. It is 
thought that A0620-00 has an outburst recurrence time of 
58 years, and more than half of that time has passed since 
1975, so the system should be gearing up for a new outburst. 
In the disc instability model of outbursts, the disc steadily 
recharges between outbursts, growing in size and density 
until it becomes unstable. In the context of this model, it 
is quite possible that since 1994, we have actually watched 
the disc expand towards the 3:1 resonance (r=0.667?Li) and 
become distorted. If this resonance is inside the disc, then 
the distortion tim e-scale is g~^Porb, or about ninety days 
(|Frank et al.ll200^ ) . If not, the process is slower, but a change 
over twelve years for A0620-00 seems reasonable. 

Now if the disc is actually circular, the Doppler map 
should be radially symmetric about the point (0,--fLi), and 
by locating the centre of th e disc, we can ident i fy the radial 
velocity of the black hole (jSteeghs fc CasaresI [20o3 ) . Oth- 
erwise, the eccentricity of the disc should be evident in a 
discrepancy between the observed and predicted locations 
of this point. We therefore implemented a search for the 
disc centre of symmetry, starting at the predicted point and 
extending ±200 km s~^ in 14 and Vy, subtracting the sym- 
metric part, squaring, and computing the mean and stan- 
dard deviation. The point with the lowest mean corresponds 
to the centre of symmetry. We iterated our search, updating 
the centre and improving the resolution until we found the 
minimum. 

Since our disc is very structured, it was difficult to find 
a region of the image unaffected by bright spots, so we per- 
formed our search on a smoothed version of the map. We 
then ran 10,000 Monte Carlo simulations to find the centre 
of symmetry, using the standard deviation of the residuals 
as the uncertainty to be sampled. Combining the results for 
our Ha and H/3 maps, we find the center of symmetry at 
(80,-220)±(40,20) km s"\ weU outside the uncertainty in 
Ki. Since we cannot explain a factor of nine of implied in- 
crease in q, we take this result, coupled with the bright spot 
locationa and the crescent features in the maps, as evidence 
for the non-zero eccentricity of the disc. 

Finally, we recognize that as the disc precesses, its 
apparent centre of symmetry should move, as this point 
roughly corresponds to the mean radial velocity of the disc. 
We made Doppler images of the first night and the two 
halves of the second night, and performed the above search 
for the centre of symmetry for each. However, the low phase 
coverage of these maps results in large uncertainties, and we 
are unable to determine a trend. The feasibility of such a 
measurement also depends on the portion of the disc par- 
ticipating in eccentricity and precession, and the exact pre- 
cession period, which is generally estimated at ~ 30 orbits. 
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Figure 6. Ha Doppler maps. Notice in the middle left panel that the bright spot is located between the gas stream trajectory and the 
Keplerian velocity of the disc along the stream, and the brighter crescent-shaped features. 
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Figure 7. H/3 Dopplcr maps. The bright spot is again between the gas stream and the Keplerian disc velocity, and the crescents are 
obvious. The gas stream itself is clearly visible in the asymmetric part (middle right). 
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Figure 8. H7 Doppler maps. Only the bright spot is prevalent. 
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Figure 9. He I 6678 Dopplcr maps. Again, only the bright spot is apparent. 
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With improved spectral resolution, a sequence of 3-4 full 
nights (for sufficient phase coverage per night) sampling the 
precession period should allow the motion of the center of 
symmetry of the disc to be resolved. 



4.3 Modulation tomography 

It is also possible to relax the assumption of Doppler tomog- 
raphy that the source flux is constant throughout the orbit 
ijSteeghslliooi ). The new technique of modulation tomog- 
raphy allows not only the imaging of average line emission 
from an accretion disc, but also maps harmonic variations on 
the orbital period. The technique is robust and flexible, but 
requires somewhat better S/N than standard tomographic 
imaging. Therefore we only consider our Ha profile here. 
The process of creating modulation maps is quite similar t o 
standard tomography, and is described in ISteegt3 (|2003h . 
The map is presented in Figure 10. 

It is clear from the fitted data (top right), that modula- 
tion tomography does a much better job reproducing the 
trail than standard Doppler mapping. Whereas standard 
mapping could only reach = 35, we were able to attain 
= 9.75 with modulation tomography. The poor is due 
to flickering. The constant portion of the map (middle row, 
left panel) looks quite similar to the maps presented in previ- 
ous sections, though the bright spot is significantly smaller. 
We still place it off the gas stream trajectory, but it is much 
less diffuse, and is located, as before, near r — 0.6i?Li. It 
is clear from this image that the inner edge of the disc in 
velocity space, and thus the outer edge of the disc in phys- 
ical space, extends as far as 0.7i?Li. The crescents are not 
nearly as prevalent in the constant part, but are apparent 
in the modulation maps. The crescent placed at 7 o'clock 
in the standard map appears in the costp map (bottom left) 
and the crescent placed at 2 o'clock appears in the sin tp map 
(bottom right). The full modulating amplitude is shown in 
the middle row, right panel. 

For an erratic source like A0620-00, these maps must 
be interpreted carefully. Modulation tomography, as men- 
tioned, maps only harmonic variations, and the demon- 
strated flickering is hardly harmonic. However, we have 
shown the eccentric disc to be a viable explanation for the 
observed phenomena, and the regions responsible for the 
crescent emission are not constant by any means. It is clear 
from the constant image that the inner edge of the disc in 
velocity space, and thus the outer edge of the disc in physi- 
cal space, extends as far as 0.7-Rli, lending credence to the 
3:1 resonance argument. For typical disc precession, super- 
humps are observed with periods 3 — 5 per cent off the orbital 
period (|Whitehurst fc Kina,199lb . If indeed this modulating 
emission is caused by viscous dissipation at the outer edge 
of an eccentric precessing disc, there should be some power 
at Porb, though a map of variation at some superhump pe- 
riod Psh would reveal larger amplitudes. As yet, such a map 
is not possible. None the less, modulation tomography has 
provided yet another piece of evidence for our eccentric pre- 
cessing disc model, and the fact that there is a component 
of disc light modulating at the orbital period could explain 
the observed light fraction from the secondary. 



5 DISCUSSION 

We have presented spectroscopic analysis of the black hole 
binary A0620-00. We measure an absorption-line radial ve- 
locity K2 = 435.4 ± 0.5 km s^^. With two measurements of 
the rotational broadening of the secondary, we flnd a mass 
ratio of g = 0.060±0.004 and a minimum mass of 3.10±0.04 
Mq for the primary obje c t. Wi th the most likely inclination 
of 41° from lGelino et al.l (|200ll ). measured in J, H, and K, 
the black hole has a mass of 11.1 Mq. The strong infrared 
flickering discussed earlier, in conjunction with unexplained 
smooth variability in the lightcurve and uncertainty in the 
disc spectrum itself, makes it difficult to estimate the true 
uncertainty in the inclination. The range of reported inclina- 
tions, 31° to 70.5°(Gelino et al. 2001 and references therein), 
results in a black hole mass between 3.7 Mq and 22.7 Mq. 
Until the nature and variability of the light from the disc 
is revealed in full, this conservative error estimate must be 
sufficient. 

We also find that the secondary contributes 44 ± 7 per 
cent of the light near 5500 A. As this means that the disc 
contributes a significant fraction of the light, especially in 
emission line regions, it becomes important to assess the 
variability of the disc, particularly if the inclination is to be 
determined by lightcurve modelling. As noted, it is common 
to assume that the disc is a constant source of light. While 
it may be valid when / is large, three observational points 
cast doubt on this assumption: 

(i) S04 performed a detailed study of flares from A0620- 
00, which in their observations have amplitudes nearing 20 
per cent of the source flux. 

(ii) Measurements of the fraction / of light contributed by 
the secondary have not been consistent. MR86 found 40± 10 
per cent at 5100 A, MR W94 found 94 ± 3 per cent at Ha, 
and lGelino et al.l (|200j ) found / > 97 per cent in J, H, and 
K, assuming a constant diluting source of light. 

(iii) Observations in all four Spitzer bands, taken approx- 
imately two weeks before our observations on the Clay tele- 
scope, show strong flickering which is highly correlated with 
simultaneous R-band lightcurves from the 1.2m telescope on 
Mt. Hopkins. 

While McClintock, Home, and Remillard (1995) rightly 
point out that the absorption line strength of the template 
will affect the observed dilution fraction, unless it can be 
shown that the measured fraction is strongly correlated with 
template line strength, a physical origin for this variation 
cannot be ruled out. Future studies could assess the true 
dependence of / on template star, as well as the long-term 
variability of /, by observing both BS 753 (MWR94's tem- 
plate) and HD 7142. Our measurements of the Ha equivalent 
widths, when compared to those of MRW94, suggest a physi- 
cally real origin for the variation, because equivalent widths 
are independent of the template star and the instrument. 
So there appears to be evidence for a disc which is brighter 
relative to the secondary than it used to be. 

Indeed, the modulations of the equivalent width even 
point to a physically non-uniform flickering, because the disc 
emission lines vary relative to the non-stellar continuum. We 
can tentatively identify the flickering with the crescents, nei- 
ther of which was present in 1994, so our conclusion seems 
viable. The line emission and the continuum may have dif- 
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Figure 10. Ha modulation map. The observed data (top left) are well-reproduced by the fitted data (top right). The middle left panel 
shows the constant part of the disc, and the middle right panel shows the amplitude of modulation. The bottom row, left and right, are 
the cosine and sine components of variation, respectively. The crescent regions arc obviously modulated at the orbital period. 
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ferent radial emissivity dependencies, which could result in 
slower modulations. We also detect several bright regions 
in the disc: one near the gas stream impact point, and two 
crescent-shaped regions on opposite sides of the disc. The 
reality of these features, as well as the non-uniform flick- 
ering, are confirmed by modulation tomography of the Ha 
disc line, which reveals variation near the bright crescents. 

First noticed by S04, the crescent-like features may in- 
dicate an eccentric disc, which is predicted for systems like 
A0620-00 with small mass ratios and large discs. It seems 
that we have observed an eccentric disc, but let us consider 
the evidence. From our observations, the following are clear: 

(i) The disc is bright and variable. The brightness is ev- 
ident in the increased dilution of the secondary spectrum, 
and the variability is clear from a number of phenomena. 
First, the trailed Ha line shows clear evidence of flickering 
events. Second, the subtracted equivalent width of the same 
line is variable beyond explanation by noise alone. Third, 
the phase-resolved light fraction cannot be reproduced by 
ellipsoidal variability on top of a constant source of light. 

(ii) The disc extends to the 3:1 tidal resonance. This is a 
simple point, clear from the Doppler and modulation maps. 

(iii) The disc is not centered on the radial velocity of the 
black hole. 

(iv) The disc is not radially symmetric, but characterized 
by bright crescent-shaped regions. 

(v) The crescent regions are modulated at the orbital pe- 
riod. 

Each of these points alone would be insufficient evidence 
to conclude that the accretion disc is eccentric and precess- 
ing. But with the exception of a direct image of the ellipti- 
cal disc, we can present a complete and coherent argument 
that this is the case. The disc has grown to tidal resonance, 
where the enhanced disc viscosity results in bright and vari- 
able rims of extra dissipation. We then observe crescents 
of extra dissipation at relatively low velocities, as expected. 
Given the viscous effects, it is predicted that the disc will 
receive a gravitational torque from the secondary, and begin 
to precess. The asymmetries introduced here shift the ve- 
locity center of disc emission away from the black hole, and 
we find that the disc is not centered on the black hole in 
velocity space. Furthermore, given the beat period between 
the precession and orbital motion, the regions of viscous dis- 
sipation should be modulated at roughly the orbital period. 
Modulation tomography reveals this to be the case. 

In retrospect, knowing that portions of the disc are 
modulated on the orbital period, we look closer at the frac- 
tion of light contributed by the secondary star, and see that 
it is not well fit by ellipsoidal modulations for a system at 
the inclination of A0620. But if another component of the 
system was variable on the orbital period, as we have ob- 
served the disc to be, then there is no need for concern. The 
physical picture, a precessing elliptical disc torqued by the 
secondary star, predicts and produces all the phenomena we 
have discussed in our data, which are of high quality. 

To put it another way, the eccentric disc hypothesis is 
nicely self-consistent. It explains why and how the accretion 
disc has changed, allows the disc to be large enough for the 
growth of eccentric modes, and predicts the phenomena that 
we observe. It is unfortunately not possible at this point 
to make an estimate of the disc eccentricity. ISmith et al.l 



(|2007l ) have shed a great deal of light on the evolution of 
disc eccentricity and energy dissipation with 3D SPH sim- 
ulations. They find that systems with q between 0.08 and 
0.24 develop low-mass eccentric discs withsuperhumps; for 
q = 0.0526, the disc exhibits a short-lived superhump and 
decaying eccentricity. All mass ratios show enhanced dissi- 
pation in the disc from the thermal-tidal instability, even 
without the eccentric modes. 

Since we have not observed a superhump, we cannot 
place A0620-00 in either category. If it falls in the more 
extreme group, the disc eccentri city is likely zero (reached 
after about 300 orbital periods) (|Smith et al.ll2007l '). In that 
case, the steady state is a massive disc. If the steady state 
is very long-lived, and the disc continues to grow, this could 
explain the enormous intensity of novae like A0620-00. If it 
fits among the less extreme mass ratios, the disc eccentric- 
ity is around 0.1-0.2, and a superhump should be observable 
with better photometry and a longer baseline (|Smith et al.l 
120071 ). A0620-00 may also be at a transition between those 
cases, and its evolution might be somewhat more erratic, 
as suggested by the SMARTS data discussed earlier. For 
example, it may toggle between states of quiescence, super- 
humps, and variability (like what we have observed here). It 
might, then, be erroneous to interpret this recent increase 
in brightness as the build towards outburst. 

While we have strong evidence that the accretion disc 
around the black hole has grown out to the tidal distortion 
radius, evolved into an eccentric disc, and started to precess, 
further study is required to verify our conclusion. Data from 
SMARTS, FLWO, and Spitzer will further quantify fiicker- 
ing, and may reveal a superhump, or some new period con- 
sistent with our results, and future programs of tomography 
will track the evolution of the accretion disc. In anticipa- 
tion of the impending outburst, and in light of progress in 
simulations, we suggest that this well-studied system not be 
disregarded or ignored, for it affords us the opportunity to 
watch the evolution of an accretion disc from quiescence to 
outburst, and the chance to test models for disc instabilities 
in X-ray novae. 
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